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Abstract: In a recent paper arXiv:1107.5048, we discussed the correlation between 
the elastic neutralino-nucleon scattering cross section, constrained by dark matter direct 
detection experiments, and fine-tuning in the electroweak symmetry breaking sector of 
the Minimal Supersymmetric Standard Model (MSSM). Here, we show that the cor- 
relation persists in the Next-to-Minimal Supersymmetric Standard Model (NMSSM), 
and its variant, A-SUSY. Both models are strongly motivated by the recent discovery 
of a 125 GeV Higgs-like particle. We also discuss the implications of the recently pub- 
lished bound on the direct detection cross section from 225 live days of XENON100 
experiment. In both the MSSM and the NMSSM, most of the parameter space with 
fine-tuning less than 10% is inconsistent with the XENON100 bound. In A-SUSY, on 
the other hand, large regions of completely natural electroweak symmetry breaking are 
still allowed, primarily due to a parametric suppression of fine-tuning with large A. The 
upcoming XENON IT experiment will be able to probe most of the parameter space 
with less than 1% fine-tuning in all three models. 
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1 Introduction 



Several experiments around the world are currently attempting to observe dark matter 
via "direct detection", i.e. measuring recoils of detector nuclei following their collisions 
with ambient dark matter particles. While no convincing observation has been reported 
so far, the sensitivity of the experiments is rapidly increasing. Currently, the best upper 
bounds on the cross section of elastic, spin-independent dark matter-nucleon scattering 
in the 10 GeV-TeV mass range come from the XENON100 experiment [1, 2], and 
are of the order 10 -44 — 10 -45 cm 2 . This is the range where the predictions of many 
attractive theoretical models of dark matter lie [3]. The most studied of these is Im- 
parity conserving supersymmetry, where the lightest neutralino Xi generically has the 
right properties to explain dark matter. It is therefore important to understand the 
implications of the direct detection bounds for supersymmetric dark matter. 

In a recent paper [4] , we pointed out a strong correlation between direct detection 
cross sections and naturalness of electroweak symmetry breaking (EWSB) in the Min- 
imal Supersymmetric Standard Model (MSSM) 1 : MSSM parameter points with lower 
direct detection cross section have more finely-tuned EWSB. This conclusion seems very 
general: It does not depend on the details of the SUSY-breaking mechanism and high- 
scale physics, and for most of the parameter space (except pure-Higgsino dark matter) 
it does not require imposing the thermal relic density constraint, so that it applies even 
in models with non-standard early cosmological history. The only assumption is the 
absence of accidental cancellations among physically distinct contributions to the dark 
matter-nucleon scattering amplitudes (e.g. Higgs- and squark-exchange diagrams). 

The main goal of the present paper is to extend this study to the Next-to-Minimal 
Supersymmetric Standard Model (NMSSM). (For reviews, see, for example, Refs. [7, 
8].). Our main motivation to study the NMSSM is the recent discovery of a new 
particle, with a mass of about 125 GeV and properties consistent with the Standard 
Model (SM) Higgs, at the LHC [9]. In the MSSM, a Higgs particle of this mass can only 
be accommodated at a price of severe fine-tuning in the EWSB [10]. Fine-tuning is 
significantly alleviated if an extra singlet (with respect to SM gauge groups) superfield 
is coupled to the Higgs sector, as in the NMSSM [10]. In this respect, a particularly 
promising variation of the NMSSM is the "A-SUSY" [11], characterized by a large 
superpotential coupling A between the singlet and doublet Higgs fields. The addition 
of a singlet has a non-trivial effect on dark matter phenomenology due to the possible 
admixture of the singlino in Xi-> as we ^ as the additional SM-singlet Higgs state. Does 
the correlation between direct-detection rates and EWSB fine-tuning persist in the 
NMSSM and A-SUSY? 

1 Related discussions have also appeared in Refs. [5, 6]. 
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(a) (b) 

Figure 1. Feynman diagrams contributing to spin-independent elastic scattering of neu- 
tralino dark matter off a nucleon in the MSSM. 



This paper is organized as follows. We begin by reviewing and updating the MSSM 
results of Ref. [4] in Section 2, including the latest results from XENON100 [2]. We then 
briefly review the structure of the NMSSM and A-SUSY in Section 3. Section 4 provides 
definitions of fine-tuning, and discusses a suppression of fine-tuning at large values of 
A, which will be important for understanding our results. The analysis procedure is 
outlined in Section 5, and the results are presented and discussed in Section 6. Finally, 
we conclude in Section 7 with a brief summary of the main results and directions for 
future work. 

2 Naturalness and Direct Detection in the MSSM: a Mini- 
Review and Update 

Spin-independent scattering of the neutralino off a quark in the MSSM is mediated 
by i-channel Higgs or s-channel squark exchanges, Fig. 1. As explained in Ref. [4], 
any cancellation between the two contributions should be regarded as accidental. We 
assume that no such accidental cancellations occur, so that the total cross section is at 
least of the same order as the contribution of either diagram class by itself. Since we 
are interested in the lowest direct detection cross section possible for a fixed amount 
of fine-tuning, it is sufficient for our purposes to examine a single diagram class. We 
focus on the Higgs exchange diagrams, since they contain the parameters that directly 
affect electroweak symmetry breaking, and we ignore the squark exchanges. The cross 
section then depends on just five MSSM parameters: /x, Mi, M 2 , tan/3, and tua- The 
same parameters suffice to determine the degree of fine-tuning in the EWSB at tree 
level, as well as the composition of the lightest supersymmetric particle (LSP). This 
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allows us to study correlations among these quantities without any assumptions about 
the nature of the SUSY breaking sector, high-scale physics such as unification, etc. 
Such a study was performed in Ref . [4] , and the reader interested in the details of the 
analysis is referred to that paper. Here, we summarize the conclusions and update the 
main plots of Ref. [4] to include the newly released XENON100 cross section bound 
from 225 days of running in 2010-2012 [2]. We also use the updated values of nuclear 
form factors, which enter the direct detection cross section predictions (from [12]): 

= 0.023, /g = 0.033, /ff = 0.26. (2.1) 

It should be kept in mind that there is significant uncertainty in the values of the form 
factors, especially that of the strange quark, which plays a dominant role for a typical 
MSSM point. 2 In particular, lattice QCD results [14] suggest a lower value ~ 0.05. 
Correlations among physical quantities (e.g. direct detection cross section and fine- 
tuning) are independent of the particular value of ff^. However, it does play a role in 
the interpretation of cross section bounds, introducing an order-one uncertainty. We 
will return to this issue in Section 6. 

The main results of Ref. [4] are summarized by the two plots in Fig. 2. First, 
there is a correlation between the direct detection cross section and the "purity" of the 
LSP, defined as p = min(F#, 1 — where F H is the Higgsino fraction. A typical 
mixed gaugino/ Higgsino LSP has a direct detection cross section around 10 -44 cm 2 ; 
significantly lower cross sections are only possible for pure gaugino or pure Higgsino 
states. There is a direct correlation between the cross section and purity of the LSP: 
Lower cross sections require higher purity. The left panel of Fig. 2 shows that almost 
no points with purity above 0.2 are allowed by the 2012 XENON100 bound. Second, 
in the case of predominantly gaugino LSP, there is a correlation between purity and 
fine-tuning of EWSB: suppressing the Higgsino admixture requires raising /i, which 
in turn increases fine-tuning. (A more detailed discussion of fine-tuning will be given 
in Section 4.) Thus, in the gaugino LSP region, points with direct detection cross 
sections significantly below 10 -44 cm 2 are necessarily fine-tuned, and the smaller the 
cross section, the more severe the required fine-tuning is. The right panel of Fig. 2 shows 
this correlation (points with greater than 90% higgsino fraction have been removed), 
and the implications of XENON100. Very few points with fine-tuning better than 1/10 
survive the 2012 XENON100 constraint, and the majority of those have LSP masses 
below 50 GeV. XENON-IT will be able to probe fine-tuning down to 1% level for most 
of the relevant mass range. 

2 We refer to [6] for a more extended discussion of the effects of varying the strange quark form 
factor. A recent discussion of the determination of this form factor can be found in Ref. [13]. 



-4- 




kJlsp [GeV] '«lsp |GeV| 



Figure 2. Left panel: Direct detection cross section vs. the LSP mass, for MSSM points 
with purity above 0.2 (red), between 0.1 and 0.2 (orange), 0.01 and 0.1 (green), and 10 -3 
and 0.01 (cyan). Right panel: Direct detection cross section vs. the LSP mass, for MSSM 
points with gaugino-like LSP. Red, green and cyan points correspond to EWSB fine-tuning 
in the intervals (0,10); [10,100); [100,1000), respectively. The black/solid and blue/dotted 
lines correspond to the XENON100 100 days/2011 [1] and 225 days/2012 [2] exclusion limits, 
respectively. The red/dashed line shows the projected sensitivity of XENON1T [15]. Real 
values of the scanned MSSM parameters are assumed, and points with strong accidental 
cancellations have been discarded. 

All above statements apply without assuming that the relic density of neutralinos is 
set by thermal decoupling, and so are remarkably insensitive to cosmological evolution 
assumptions. In the Higgsino region, there is no correlation between direct detection 
cross section and fine-tuning. However, if one further assumes standard cosmology 
and thermal decoupling, the correct relic density for Higgsinos requires /i ~ 1 TeV, 
corresponding to EWSB fine-tuning of about 1/500, independent of the direct detection 
bounds. 

3 The NMSSM and A-SUSY 

The NMSSM is obtained by adding to the MSSM a gauge-singlet superfield 5, with a 
superpotential containing 

XSH U -H d + ^S 3 . (3.1) 

This is the most general superpotential consistent with a Z 3 symmetry under which 
5, H u and fields have charge 1/3. The fiH u Hd term present in the MSSM is not 
allowed by this symmetry. The soft SUSY-breaking Lagrangian contains the following 
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new terms: 

m 2 s \S\ 2 + (XA X H U • H d S + -kA k S 3 + h.c.) , (3.2) 

while the B/LiH u Hd term of the MSSM is forbidden by the Z 3 symmetry. When the 
singlet field receives a vacuum expectation value (vev), s = (S), after SUSY breaking, 
effective \i and B terms are generated: 

/i e// = As, B eff = A x + ks . (3.3) 

The two bosonic components of S give one CP-even and one CP-odd Higgs bosons, 
which mix with the MSSM Higgses to produce 3 CP-even and 2 CP-odd mass eigen- 
states. The extended Higgs sector is described in terms of the following parameters: 

A, k, A X) A K , m 2 Hu , m 2 Hd , m 2 s . (3.4) 

Of these, A is particularly important in accommodating a 125 GeV Higgs, since it enters 
the tree-level mass of the CP-even Higgses. Large values of A are needed to raise the 
Higgs mass. Renormalization group evolution drives A to larger values at high energies, 
and eventually A hits a Landau pole. Requiring this pole to lie above the conventional 
grand unification (GUT) scale leads to an upper bound on the weak-scale value of A, 

A < 0.75. (3.5) 

In this paper, we refer to the theory with A obeying this condition as the "NMSSM". 
However, it is also interesting to consider the regime of larger A. Even though the 
simple prediction of gauge coupling unification is lost, such theories provide perfectly 
consistent descriptions of physics at the weak scale, and require the least fine-tuning 
among models of this class to incorporate the 125 GeV Higgs [10]. Avoiding a Lan- 
dau pole below 10 TeV (which would almost certainly be inconsistent with precision 
electroweak constraints) requires 

0.75 < A < 2.0. (3.6) 

In this study we will refer to the theory with A in this range as the "A-SUSY" [11]. 
Above the 10 TeV scale, A-SUSY needs to be incorporated into a more fundamental 
ultraviolet (UV) theory; see, for example, Refs. [16, 17] for recent attempts at UV 
model building in this context. 

The fermionic component of S ) the "singlino", mixes with the four neutralinos 
of the MSSM to give 5 neutralino mass eigenstates Xi- The lightest of these, Xi> 18 
assumed to be the LSP and is stable due to conserved R-parity. We will assume that 
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all of dark matter present in the Earth's neighborhood consists of the Xi particles. The 
neutralino mass matrix has the form 

( Mi —m z s w cos (5 m z s w sin (3 \ 

M 2 rnzCw cos f3 —m z c w sin /3 

—m z s w cos /3 mzCw cos /3 ~l Jj eff —\vs'mf3 

m^s^sin/3 — m^c^ sin /3 ~f 1 eff — Xvco8/3 

\ — A'usin/? — A^cos/3 ^j^eff J 

(3.7) 

where and are the sine and cosine of the Weinberg angle, and tan/3 = v u jv^. 
The only two new parameters that enter, beyond those listed in (3.4), are the weak-ino 
soft masses M\ and M 2 . 

3.1 Direct Detection 

The spin-independent LSP-nucleon scattering in the NMSSM and A-SUSY occurs via 
the same diagrams as in the MSSM, see Fig. 1, except now all three CP-even Higgses 
can be exchanged in the i-channel. We again ignore the squark diagram, with the same 
motivation as in the MSSM study. The scattering amplitude can be easily computed 
using the NMSSM Feynman rules listed, for example, in Ref. [8]; in particular, the 
LSP-Higgs coupling is given by 

9(hiXiXi) = -^ s H d nH u n~ s + s Hu n s n~ s + ssn^n^) - s s n~ s n~ s 

+ j( s H d riBnH d - s Hu n § n Au ) - ^{s Hd n w n &d - s Hu n^n &u ) , (3.8) 

where s a and n# denote the relevant components of the z-th Higgs mass eigenstate 
(i — 1 . . . 3) and the lightest neutralino, respectively. The cross section is then obtained 
using the standard formalism [3]; we use the values of nuclear form factors listed in (2.1). 

4 Fine-Tuning in MSSM, NMSSM and A-SUSY 

The supersymmetric parameters must produce the weak scale; this is the origin of 
potential fine-tuning, reflected in the tree-level relation for the Z boson mass: 

1 9 m ?f, — tan 2 /3 m 2 M 9 

This relation holds in the MSSM as well as in the NMSSM and A-SUSY. A reasonable 
expectation is that /i and rriH u must also be around the mz scale, and some fine-tuning is 
required to obtain the correct mz if these terms are far above this scale. A quantitative 
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A max Aj , A^ = 



(4.2) 



measure of the amount of fine-tuning is the variation in m 2 z resulting from variations 
of the fundamental Lagrangian parameters of the theory. This can be defined as 

d log m| 
d log & 

where the index i runs over all independent parameters in the Lagrangian. There is an 
ambiguity over the scale at which the parameters & are defined; a common choice in 
the literature is the GUT scale [18]. We will instead use ^ defined at the weak scale, 
since we are working with general (N)MSSM without assuming specific SUSY-breaking 
scenarios. 3 At tree-level, the only relevant parameters in the MSSM are /x, m^, rriH d 
and b] analytic expressions for A^, i = 1 ... 4 are given in Refs. [4, 19]. 

In the NMSSM, minimization of the Higgs potential relates the scalar vevs v U) Vd 
and s to the Lagrangian parameters [7, 20]: 

2 

Ei = m 2 Hu +fi 2 + 2X 2 vj + 9 -{v 2 u - vj) - V ^^{A X + ks) = , 

2 

E 2 = m 2 Hd + /x 2 + 2\ 2 v 2 u + 9 -(v 2 - vl) - —fi(A\ + ks) = , 

2 v d 

E 3 = m 2 + ^A Kf i + 2^/i 2 + 2A 2 (^ 2 + v 2 ) - 2\KV u v d - \ 2 v u v d — = , (4.3) 
A A z fi 

where 

m 2 z = g 2 v 2 , v = ^vl + v 2 d = 174 GeV, (jl = \s. (4.4) 

We defined g 2 = (g 2 + g 2 )/2 « 0.52, where ^ and # 2 are the SM C/(l) y and SU(2) L 
couplings respectively. We ignored loop corrections, but it should be kept in mind that 
these can be important; the first equation, in particular, receives a large correction 
from top and stop loops. 

The variations of rrtz under input parameter changes can be calculated in the 
following way [20]: Imposing that the minimization conditions (4.3) continue to hold 
under variations of the input parameters, one obtains 

= E + + a^*"/> + f * = • < 4 ' 5 ' 

for i — 1 ... 3, where j runs over the fundamental parameters listed in Eq. (3.4). These 
three equations can be solved for 8m 2 Z) 5tan/3, and 5fi. Defining 

dE j „ dEi „ dEi „ dEi „ tA N 

wr p "< 8^r z - adrp =T - ^ =M " <46) 



3 It is possible that a cancellation that appears finely tuned from the weak-scale point of view may 
be rendered natural in a particular high-scale theory; this would be missed by our approach. 
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we obtain 



d log m| 




6 5m| 




d log & 









6 E iH ^ W 



(4.7) 



The fine-tuning is then calculated by taking the maximum of these A$'s, as described 
earlier. 



4.1 Fine- Tuning Suppression in A-SUSY 

Since we define fine-tuning as sensitivity with respect to weak-scale input parameters, 
no renormalization group evolution is involved in computing it, and the formalism 
described above for the NMSSM applies equally well to A-SUSY. However, there is a 
very interesting, and potentially important, parametric suppression of fine-tuning at 
fi 3> vrtz and mn u 3> mz in A-SUSY regime that is not present in the NMSSM. We 
will discuss this phenomenon in this subsection. 

In the MSSM, the sensitivity of vrtz to /i is given simply by 



4^ 2 



m 



2 ' 



(4- 



so that fine-tuning scales as (/i/ra^) 2 at large /i, independent of all other parameters. 
Likewise, large values of 7Uh u give fine-tuning of order (ra^/m^) 2 . In the NMSSM and 
A-SUSY, the situation is more complicated, since /i is not an input parameter, but an 
output of the minimization of the scalar potential. In particular, cancellations between 
terms in Eq. (4.1) may occur naturally ) if the minimization conditions force a particular 
relation of \i to 7Rh u and/or rriH d - This is exactly what happens, quite generically, in 
the regime A>g: the fine-tuning at large \i and mn u is parametrically suppressed by 
a factor of (g/X) 2 . This is extremely relevant for the LHC searches for supersymmetry: 
Large A allows one to raise rriH u without increased fine-tuning; the stop masses, which 
control the size of the quantum correction to ra# u , can then be raised by roughly the 
same amount. If X/g ~ 3, as is possible in A-SUSY, stop masses in excess of 1 TeV can 
be completely natural [10]. 

To understand the parametric suppression of fine-tuning at large A, consider the 
model in the particularly simple limit, A\ = A n = 0. In this limit, the minimization 
conditions (4.3) can be solved analytically (see Appendix A for details). The result is 



/i 



rriy = 



m 



rric 



m 



m 



h 



777; 



mi 



m 



m 



(4.9) 
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Figure 3. Fine-tuning vs in the NMSSM (left) and A-SUSY (right). 

where corrections of higher orders in g 2 /A 2 have been ignored, and the functions f\ 
and f% are given in Appendix A. These functions are parametrically 0{1) if all their 
arguments are order-one numbers, so that generically, 

2 2 

2 # 2 # 2 r a -x n\ 

m z ~ ^m Hu ~ — /i , (4.10) 

making the large- A suppression of fine-tuning manifest. 

For general A\ and no analytic solution can be found, but there is still a simple 
way to argue that the suppression of fine-tuning at large A should persist. Consider A 
dependence of the parameter sensitivities A^ in Eq. (4.7). While Ti and Mi appear in 
both numerator and denominator, appears only in the numerator and Zi only in 
the denominator, so the magnitude of A^ is set roughly by the size of ^P^-'s relative to 
the m|Zj's. All Zj's contain a factor A 2 /g 2 ) a term not found in the i^/s, hence A^'s 
are suppressed by (g/X) 2 . 

To further illustrate this phenomenon, we plot \i vs. the amount of fine-tuning in 
the NMSSM and A-SUSY in Figure 3. The coupling A at each point varies randomly 
between 0.4 and 0.75 in the NMSSM scan, and between 0.75 and 2.0 in A-SUSY scan. 
Other details of the scans (which will also be used to study direct detection/fine-tuning 
correlation) are described in Section 5. Side- by-side comparison of the plots clearly 
shows that lower fine-tuning for the same value of ji is possible in A-SUSY, consistent 
with the arguments given above. 

5 Analysis Setup 

To study correlations between direct detection cross section, fine-tuning and other 
physical quantities, we performed a scan over the NMSSM/ A-SUSY parameters. The 
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fine-tuning and direct detection cross section (ignoring the squark contribution and all 
loop corrections) are completely determined by the 7 Higgs-sector parameters listed 
in Eq. (3.4) and the weak-ino soft masses Mi and M 2 . For convenience, we used 
Eqs. (4.3) to interchange Tn 2 Hu ^m 2 Hd ^ and m 2 s with tan/3, /i, and m%. With mz fixed, 
there are 8 parameters to be scanned over. The boundaries of the scan are as follows 
(all dimensionful parameters are in GeV and can have either sign): 

• fi: (70, 10000). 

Lower values of /i result in light charginos, failing the LEP-2 bound [21]. Points 
with /i above 10 TeV are fine-tuned at the 10 -3 level or worse. 

• A: (0.4, 0.75) in the NMSSM scan, (0.75, 2.0) in A-SUSY scan. 

The lower bound is chosen to get a sizable F-term contribution to the tree level 
Higgs mass, necessary to obtain a 125 GeV Higgs without significant fine-tuning. 

• k: ±(10- 5 ,0.65). 

The upper bound is the perturbativity bound in the NMSSM (i.e., k remains 
perturbative up to the GUT scale). Since the singlino mass term is ms = 2/^/x/A, 
the region k < 10 -5 typically results in very light (<10GeV) singlino LSPs, 
which are not constrained by XENON100 and will not be considered here. 

• tan/3: (1.7, 54). 

The upper limit is set by the perturbativity bound on the bottom Yukawa cou- 
pling [22]; in addition, the F-term contribution to the Higgs mass is suppressed 
at large tan/3, so that region is disfavored the 125 GeV Higgs. 

• A X ,A K : (1, 10000). 

• Mi. (10, 10000), M 2 : (80, 10000). 

We generate points randomly distributed with uniform weight in A, ft, and tan/3, 
and randomly distributed with logarithmic weight in the other five (dimensionful) pa- 
rameters. We impose the following constraints: 

• The lightest chargino must be heavier than 103 GeV to satisfy the LEP-2 bound [21], 
and must be heavier than the lightest neutralino. 

• The LSP mass must be above 10 GeV (since the current XENON100 data is not 
sensitive to lower LSP masses). 

• No tachyonic Higgses in the CP even or CP odd sector. (All Higgs masses are 
computed at tree level.) 
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• There is at least one CP-even Higgs with a tree-level mass between 100 and 150 
GeV and a singlet fraction of at most 50%. These requirements select the part 
of parameter space roughly consistent with the 125 GeV Higgs discovery. 

• The LSP contribution to the invisible width of the Z must be less than one 
standard deviation of the measured neutrino contribution: T^ xx <4.2MeV if 
m x <mz/ i 2 [21, 23]. The light LSP candidates that survive this constraint are 
mostly bino. 

After imposing these conditions, about 30000 points each remain in the NMSSM 
and A-SUSY scans. All plots in the following section are based on these data sets. 

Note that we do not demand correct thermal relic density of the neutralino for the 
points in our scan. There are two reasons for this. First, relic density depends on many 
model parameters (e.g. squark and slepton masses) in addition to the ones we scan 
over, and in most cases we expect that it can be "fixed" by an appropriate choice of 
those parameters, with no effect on EWSB fine-tuning. Second, not demanding that 
the LSP be a thermal relic gives our results broader applicability, including scenarios 
with non-thermal dark matter production, non-standard cosmological evolution, etc. 
Only in one special case (pure-Higgsino LSP) shall we consider the effect of including 
the relic density constraint. 

In discussing the correlations between the direct detection cross section and other 
quantities, it is important to distinguish parameter points where the direct detection 
amplitude is strongly suppressed for accidental reasons. A quantitative measure of an 
accidental cancellation is the dependence of the cross section on the scan parameters: 
An abnormally strong dependence indicates an accidental cancellation. We will use the 
"accidentally" , defined as 

(5.1) 

to quantify the presence of accidental cancellations at our scan points. This is analogous 
to our definition of fine-tuning. 

6 Results and Discussion 

First, consider the correlation between the direct detection cross section and the com- 
position of the LSP. This is illustrated in Figs. 4 and 5. As in the MSSM, an LSP with 
a "generic" composition (i.e. roughly equal mix of Higgsino, singlino and gaugino) 
is predicted to have a cross section of about 10 -45 cm 2 or higher. This is precisely 
the region now probed by XENON100. Points with cross sections below 10 -45 cm 2 
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d log & 
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Figure 4. Direct detection cross section vs. Higgsino fraction of the LSP, in the NMSSM 
(left) and A-SUSY (right). In the top row, all points are included; in the bottom row, points 
with accidental cancellations (Acc> 50) are discarded. 



either have accidental cancellations in the cross section, or have an LSP with the Hig- 
gsino fraction close to either 1 ("pure Higgsino" regime) or ("pure gaugino/singlino" 
regime). This can be easily understood by examining the LSP-LSP-Higgs coupling in 
Eq. (3.8). Schematically, the four terms in the coupling have the form HS ) HW, HB, 
and SS. The contributions of the first three terms are clearly suppressed in either 
pure Higgsino or pure gaugino/singlino limits. The fourth term is not suppressed in 
the pure-singlino limit, but its contribution is proportional to the coupling k and to 
the singlet-doublet Higgs mixing angle, both of which can be small. This qualitative 
behavior is not sensitive to A, and thus applies equally in the NMSSM and A-SUSY 
regimes. 

The dependence of the cross section on the singlino fraction of the LSP is quali- 
tatively similar: cross sections below the "generic" 10 -45 cm 2 level occur for singlino 
fractions close to or 1, see Fig. 5. In principle, small cross sections can also occur for 
an LSP with order-one singlino and gaugino fractions, and small Higgsino admixture. 
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Figure 5. Direct detection cross section vs. singlino fraction of the LSP, in the NMSSM 
(left) and A-SUSY (right). In the top row, all points are included; in the bottom row, points 
with accidental cancellations (Acc> 50) are discarded. 



However, such points require a near-degeneracy of the singlino and gaugino (bino or 
wino) diagonal terms in the neutralino mass matrix, Eq. (3.7). This is because there is 
no direct singlino-gaugino mixing entry in the mass matrix, so the mixing must occur 
via Higgsinos, which however are nearly decoupled in this region, and can only lead to 
large mixing if the corresponding diagonal entries are nearly degenerate. Such points 
rarely occur in the scan, and typically fail the accidentality cut, which explains their 
absence in the lower panel of Fig. 5. 

It is worth pointing out that the direct detection cross section can be suppressed 
arbitrarily if the LSP is a pure singlino, and A and n are both small. In this region, 
however, adding a singlet to the MSSM loses its strongest motivation, since small A 
suppresses the tree-level contribution to the Higgs mass, and incorporating the 125 
GeV Higgs becomes as difficult as in the MSSM. For this reason, we restricted our scan 
to A > 0.4, and the small- A region does not appear in our plots. In addition, in A — >> 0, 
k limit the singlino LSP is completely decoupled, so it is difficult to imagine a 
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Figure 6. Direct detection cross section vs. EWSB fine-tuning, in the NMSSM (left) and 
A-SUSY (right). In the top row, all points are included; in the bottom row, points with 
accidental cancellations (Acc> 50) are discarded. The lines represent approximate lower 
bounds from Eq. (6.1). 



production mechanism, thermal or non-thermal, that would give it the observed relic 
abundance. 

The most interesting correlation is between the direct detection cross section and 
EWSB fine-tuning. This is shown in Figs. 6 and 7. In these plots, we do not include 
points with Higgsino fraction in the LSP above 0.9; the predominantly Higgsino case 
will be discussed separately below. The correlation found in the MSSM is preserved 
in the NMSSM and A-SUSY: Points with lower direct detection cross sections have 
stronger fine-tuning, unless accidental cancellations occur. The physical origin of the 
correlation is the same as in the MSSM: suppressing direct detection cross section 
requires suppressing the Higgsino fraction, which can only happen (for a fixed LSP 
mass) by raising /i, but this necessarily increases fine-tuning. It does not matter for 
this argument whether the LSP is gaugino, singlino or some combination of the two. 
The only caveat is the pure-singlino case with A — >> and k —> 0, where cross section 




Figure 7. Direct detection cross section vs. LSP mass, in the NMSSM (left) and A-SUSY 
(right). Red, green, cyan and yellow points correspond to EWSB fine-tuning in the intervals 
(0,10); [10,100); [100,1000), and > 1000, respectively. Points with accidental cancellations 
(Acc> 50) are discarded. Lines denote XENON bounds and projections, as in Fig. 2. 



can be suppressed with no fine-tuning price; however, as explained above, this region 
is not physically interesting. 

As in the MSSM, there is an approximate lower bound on direct detection cross 
section consistent with a given amount of fine-tuning. As a rough bound, we estimate 
from Figs. 6 that, in the absence of accidental cancellations, 

10" 45 cm 2 10" 46 cm 2 
<rsi £ ^ (NMSSM); a S i £ (A - SUSY). (6.1) 

The bound depends on the LSP mass, and the above estimates apply to the lightest 
LSP masses included in our scan, ttilsp = 10 GeV; the bound is higher for heavier 
LSPs. This is illustrated in Fig. 7, which also shows current and projected XENON 
bounds. In the NMSSM, most points with fine-tuning better than 1/10 are already 
ruled out, although some remain. Most of the remaining points lie either in the very 
low LSP mass region, where the XENON bound is weakened, or in the "tail" of the low 
cross-section points in the 100 — 300 GeV mass range. The "tail" region, which has no 
counterpart in the MSSM (c/. Fig. 2), arises from the region of parameter space with 
mostly-gaugino LSP and the Higgsino and singlino masses at roughly the same scale. 
The presence of the singlino lowers the Higgsino fraction in the LSP for the same value 
of /x, resulting in lower cross section for the same level of fine-tuning. 

In contrast to the MSSM and the NMSSM, in A-SUSY a large region of com- 
pletely natural parameter space remains unconstrained by XENON100 (see right panel 
of Fig. 7). The reason for this is the parametric suppression of fine-tuning at large A, 
discussed in Section 4.1. This effect allows points with the same amount of fine-tuning 
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Figure 8. Same as Fig. 7, but with f$ = 0.05. 

to have a significantly higher \i than they would in the MSSM or the NMSSM. Higher 
H suppresses the Higgsino fraction, which in turn lowers the cross section. Overall, we 
conclude that A-SUSY allows significantly lower direct detection cross sections than 
either the MSSM or the NMSSM, without conflict with naturalness. 

When the LSP is predominantly Higgsino, the direct detection cross section in the 
MSSM can in principle be suppressed very significantly without paying a fine-tuning 
price [4]: Keeping \± ~ mz and raising Mi 3> m^, M 2 3> mz does not introduce 
tree-level tuning, and the only limitation comes from the one-loop contributions of 
weak-inos to m 2 Hu and m 2 Hd . This "Higgsino loophole" can only be removed if thermal 
decoupling is assumed, in which case \i 3> M z is required to match the observed relic 
density. In the NMSSM, the region /i ~ m^, Mi 3> m^, M 2 3> mz does not generically 
produce a pure Higgsino LSP, because the singlino does not decouple: the singlino 
mass is of order /x, and the singlino-Higgsino mixing terms are of order u, see Eq. (3.7). 
Since the mixed Higgsino-singlino LSP has a large direct detection cross section, one 
would naively expect that there is no "Higgsino loophole" in the NMSSM. However, 
although the singlino is not parametrically decoupled, a modest hierarchy between A 
and k appears to be sufficient to suppress the singlino fraction in the LSP, and hence 
the cross section, while keeping \± ~ mz- As in the MSSM, the neutralino annihilation 
cross section in this parameter region is too high to produce the observed dark matter 
abundance in the standard thermal decoupling scenario unless \i ^TeV, reintroducing 
fine-tuning. 

As already mentioned in Section 2, there is significant theoretical uncertainty on 
nuclear form factors that enter the direct detection cross sections, with the strange 
quark form factor being especially important. Throughout the paper, we used f$ = 
0.26. An alternative lattice estimate is much lower, f$ ~ 0.05 [14]. To illustrate the 
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effect of this uncertainty on our conclusions, we repeated the scan with f$ = 0.05. The 
results are shown in Fig. 8. The correlation between the fine-tuning and direct detection 
cross section is unaffected. Numerical values of cross section are scaled down by an 
(9(1) factor, so that more completely natural NMSSM points survive the XENON100 
constraint. A better understanding of the nuclear physics is clearly very important for 
clean theoretical interpretation of the direct dark matter searches. 

7 Conclusions 

In this paper, we continued our study, initiated in Ref. [4], of the implications of 
dark matter direct detection experiments on supersymmetric models. In Ref. [4] we 
found that the MSSM, for "generic" parameters, predicts a spin-independent elastic 
neutralino-nucleon scattering cross section of the order of 10 -45 — 10 -44 cm 2 (depending 
somewhat on the poorly known strange quark form factor) or higher. Cross sections in 
this range are currently being probed by XENON100, which already places meaningful 
constraints. Suppressing the cross section below the "generic" level in the MSSM 
requires (barring accidental cancellations) that the LSP be either a pure gaugino or 
a pure Higgsino. In the first case, lowering the direct detection cross section requires 
raising /i, and therefore introducing fine-tuning in the EWSB. In the second case, 
requiring that the Higgsino be a thermal relic implies EWSB fine-tuning of about 
1/500, independent of direct detection bounds. 

The recent LHC discovery of a new particle with 125 GeV mass and properties 
consistent with the SM Higgs puts significant pressure on the MSSM, since fine-tuning 
of order 0.1% is required to accommodate it in this model. This motivates considering 
supersymmetric models with non-minimal Higgs sectors, where new contributions to 
the tree-level Higgs mass can easily arise. The simplest example is the NMSSM, where 
a single gauge-singlet superfield is added, and a 125 GeV Higgs can be incorporated 
with far less fine-tuning. The reduction of tuning is especially striking in the version of 
the NMSSM with strong doublet-singlet Higgs coupling, the so-called A-SUSY. In this 
paper, we extended the analysis of Ref. [4] to the NMSSM and A-SUSY. We found that 
the qualitative correlations between the dark matter direct detection cross section, the 
LSP composition, and the EWSB fine-tuning found in the MSSM essentially persist in 
these non-minimal models as well. Numerically, the minimal cross section allowed for 
the same level of EWSB fine-tuning is somewhat decreased in the NMSSM compared 
to the MSSM, and is further decreased, rather significantly, in A-SUSY. We discussed 
the physical origin of these effects. We found that the current XENON100 cross section 
bounds are in mild tension with the MSSM and the NMSSM, excluding most points 
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with fine-tuning of 1/10 or better, while large parts of completely natural parameter 
space are still allowed in A-SUSY. 

As in Ref . [4] , we took a deliberately broad approach to the supersymmetric model 
parameter space. We do not assume any particular model of SUSY breaking; instead, 
we scan over unconstrained weak-scale parameters of each model. We do not impose 
the relic density constraints on the LSP. 4 We also do not fully utilize the recent data 
concerning the 125 GeV Higgs candidate: We only require broad consistency of the 
Higgs spectrum with the data, demanding that a tree-level mass of at least one mostly- 
doublet CP-even Higgs be between 100 and 150 GeV. Imposing any combination of 
additional constraints would select a subspace of the broad parameter space studied 
here. Stronger conclusions can be obtained with such added constraints, but they 
would be less generally applicable. It would be interesting to perform such studies in 
the future. 

Looking ahead, we can anticipate further dramatic improvement in the experimen- 
tal sensitivity of direct detection dark matter searches within a few years. Our results 
indicate that, if dark matter is supersymmetric, the searches will very likely be suc- 
cessful: Only supersymmetric models with sub-percent levels of EWSB fine-tuning, or 
accidental cancellations, will escape detection by experiments with sensitivity levels 
expected of, for example, XENON1T. In this paper, we showed that these statements 
apply not only in the MSSM but also in well-motivated non-minimal supersymmetric 
models. This underscores the importance of the continuing direct dark matter searches 
for fundamental physics. 
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A Fine- Tuning in A-SUSY: an Analytic Example 

A simple analytic solution to the NMSSM scalar potential minimization conditions, 
Eqs. (4.3), can be obtained in the absence of the two A terms, A\ = A K = 0. Define 

4 Relic density constraints in the NMSSM are well known; see, for example, the review article [7]. 
Relic density constraints in A-SUSY were considered in Ref. [24]. 
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tp = tan/3 = v u /vd- The first two equations are a linear system with respect to v 2 and 
fj 2 . Solving this system provides expressions for these quantities in terms of tp: 

2 2 1 / a2 i - 2\ 

where R = k/X^ fh 2 d = rn 2 Hd /m 2 Hu , and terms of order g 2 /A 2 were neglected since we are 
interested in the large-A limit. The last of the equations (4.3) is then used to determine 
tp: 

™ 2 s(tp- 1) + 2k 2 + m 2 d ) + l - (1 - - m\{\ - Rip 1 )) (l + t}-2Rtp) = , (A.2) 

where m| = m 2 s /m 2 Hu . This is a cubic algebraic equation on tp. For Lagrangian 
parameters consistent with EWSB, it must have a real, positive solution, such that 
ji 2 and v 2 are also real and positive, and the vacuum is stable (i.e., the Higgs mass 
matrix does not have negative eigenvalues). Eqs. (A.l) and (A.2) together yield the 
equations (4.9) in the main text, and implicitly define the functions fi and / 2 . (Explicit 
expressions can be obtained, but are not particularly illuminating.) For understanding 
fine-tuning at large A, it suffices to observe that if all coefficients in Eq. (A.2) are order- 
one numbers, the resulting tp is also generically an order-one number. An example of 
Lagrangian parameters with consistent EWSB and the parameter hierarchy of interest 
to us is 

A = 2.0, R = 0.7, fh 2 s = 2.0, fh 2 d = 0.5, m 2 Hu < . (A.3) 

At this point, tp ~ 1.63, and we obtain rn 2 z /m 2 Hu w 0.06 and m 2 z /fi 2 & 0.05, with no 
fine-tuning. (Note that this example is meant to be a simple illustration of the tuning 
suppression at large A, and is not fully realistic: in particular, it does not contain a 125 
GeV Higgs.) 
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